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Abstract 

We show that rational data of bounded input length are uniformly dis- 
tributed with respect to condition numbers of numerical analysis. We deal 
both with condition numbers of Linear Algebra and with condition num- 
bers for systems of multivariate polynomial equations. For instance, we 
show that for any w > 1 and for any n x n rational matrix M of bit length 
0(n^ log n) + logu;, the condition number k{M) satisfies k{M) < wn^^'^ with 
probability at least 1 — 2w~^. Similar estimates are shown for the condi- 
tion number ^nor-m of M. Shub and S. Smale when applied to systems of 
multivariate homogeneous polynomial equations of bounded input length. 
Finally we apply these techniques to show the probability distribution of 
the precision (number of bits of the denominator) required to write down 
approximate zeros of affine systems of multivariate polynomial equations of 
bounded input length. 
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1 Introduction 



In we exhibited examples of systems of multivariate polynomial equations 
with integer coefficients of small input length such that the required precision 
to write down an approximate zero (in the sense of is exponential in 
the number of variables (cf. also Section ^ below). These examples are 
based on the fact that the condition number finorm of [^] and HI is doubly 
exponential in the number of variables. 

In his seminal paper of 1948 (cf. |^6|), A. Turing introduced two condition 
numbers in Linear Algebra : the condition numbers M and N. A. Turing 
also showed that the M— condition number is a sharp worst case estimate 
of the loss of precision in linear Algebra Computations. The same condition 
number was also discussed by J. von Neumann and collaborators (cf. p^) 
and by J. H. Wilkinson (cf. ||, [||). 

In his paper, A. Turing already exhibited examples of systems of linear equa- 
tions with rational coefficients of small input length such that the condition 
number becomes too large. 

In all these examples, input length means the number of tape cells (also 
digits) required to represent the input in a Turing machine. 
Nevertheless, there are deep and meaningful studies on the probability dis- 
tribution of condition numbers of numerical analysis. The reader may follow 
the estimates for the Linear Algebra case in [Q , 0, |^ and the estimates 
for the Non- linear case in |3^, |^ or [^]. Most of these studies show that 
condition numbers are (almost) polynomial in the dimension of the source 
space with high probability. Hence, our examples in and Turing's exam- 
ples seem to be in contradiction with these studies. 

This "contradiction" is misleading. The probability distribution obtained in 
the series of papers [||], [0,|3|, @ is based on a "continuous model 

for numerical analysis procedures" . A "continuous model" assumes that in- 
puts belong to some "continuous source space". Consequently, studies on 
the probability distribution under a "continuous model" mean "probabil- 
ity distribution of condition numbers of inputs in some continuous source 
space" . 

For instance, in the Linear Algebra case, "inputs" are projective points in 
the "continuous source space" where Ain(M) is the vector space 

of all ra X n real matrices. Thus, if k{M) denotes the standard condition 
number of matrices M £ A4„(M) (see more details below), the "probability 
distribution" of k means the probability distribution of k when M runs over 
the real projective space P(A^„(1R))- 

On the other hand, real life computing is "discrete" . This essentially means 
that real life computers take as input matrices M G P(A^„(i^)) where K 
is some "discrete field" {K = Q, for instance). Moreover, an essential pa- 
rameter in real life computing is the input length. The input length of 
M G A^n(Q) is the number of digits required to represent the entries of 
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M in a computer (cf. also Subsection 2.4 below). Hence, real life comput- 
ing demands the probability distribution of condition numbers of numerical 
analysis for inputs of bounded length in a "discrete space" as, for instance, 
the rational projective space P(7W„(Q)). These pages are devoted to show 
this kind of estimates within a discrete context. We deal both with the 
Linear and the Non-linear cases. 

The rest of the Introduction is devoted to state the new outcomes in a 
succinct form. Subsection below shows the probability distribution of 
condition numbers k and /j, for rational matrices with bounded input length. 



Subsection 1.2 shows the probability distribution of condition number ^norm 
of M. Shub and S. Smale for systems of multivariate homogeneous polyno- 
mial equations of bounded input length whose coefficients are Gauss ratio- 



nals. Finally, in Subsection |1.3| below we show the probability distribution 
of the precision required to write down affine approximate zeros of systems 
of polynomial equations of bounded input length. 

An advised reader may realize that our statements imply that rational data 
of bounded input length are uniformly distributed with respect to condi- 
tion numbers of numerical analysis. In fact, our statements below exhibit 
concrete estimates for the discrepancy (i.e. for the difference between the 
"continuous" and the "discrete" probabilities). In other words, we show the 
uniform distribution of rational data but paying attention to A. Schonhage's 
Lemma : "Do care about constants!". 

1.1 The Linear Algebra Case 

For every square matrix M € A4„(M), let k{M) be the standard condition 
number of M. This condition number is a bound for the loss of precision of 
most linear algebra procedures. It was introduced at the end of the forties 
by several authors : A. Turing in pq|, J. von Neumann and coauthors (cf. 



[22|) and J. H. Wilkinson (cf. ||39|, |38|]). It has been extensively studied in 
[]34| , [Tsfl , 10, The condition number k{M) is given by the following 

identity 

k{M) := ||M-^||||M||, 
where ||M|| is the norm of M as linear operator, namely 

||M|| := sup{||Mu|| : u G M", ||u|| = 1}. 

The condition number k is naturally defined on the projective space 
IP(-A4n(IS))- Thus, we denote by k the mapping 

k : P(7W„(M)) — >R+. 

Let P(A4„(Q)) be the projective space defined by all n x n matrices with 
rational entries. 
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In these pages, the input length of a projective point P in P(A4„(Q)) defined 
by a matrix M £ F{M 

?i(Q)) is given as the logarithm of the Northcott— 
Schmidt height of the projective point P £ P(A^„(Q)) (of. Subsection 2.4 
below). 

The following statement holds : 

Theorem 1 (Main Theorem for Linear Algebra) With the same no- 
tations and assumptions as above, there is a universal positive constant 
ci G M, ci < 10 such that the following holds : 

Let h,w £ M be positive real numbers. Assume that w > 1 and h is big 
enough to satisfy the following inequality : 

h > cin^ log n + log w. 

Then, for every matrix M G P(A1„(Q)) of bit length at most h, the condition 
number k{M) satisfies k{M) < wn^^"^ with probability at least 

1-1. 

w 

The same assertion holds for the condition number n introduced by S. Smale 
in (cf. Subsection |2.2| below) and for the loss of precision of the gradient 
method (cf. II). 

The following constants are going to be used in what follows. Let ^ G N 
be a positive integer number. By G M we denote the volume of the 
^-dimensional unit ball in Mf. Namely, 



£r(|) 



For = 0, we write Kq := 1. Finally, let us introduce the constant S^™-) 
given by the following identity 

m-l , V 

Using E. Artin's estimates for the gamma function (cf. |Q]) the following 
estimates also hold : 

a/vT 

and 

6("^) o 

< esm 2 ei2 . 

Theorem |l] above is a direct consequence of the following more technical 
(but more complete) statement. 
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Theorem 2 With the same notations as, above let H,e gM. be two positive 
real numbers. Assume that H > 1 and e > 0. For every matrix M G 
¥{Ain{Q)) of NorthcottSchmidt's height at most H, the condition number 
k{M) satisfies k{M) > 1/e with probability at most 



B(e,n) 



where 



B{e,n) := C(n^ 



H - C{n) ' 



C(n) := C(n') 



+ 1 + 



K„2C(n2-l) 
C is Riemann's zeta function and 

Tn := (2max{4,n})(2"*+4«') 



1.2 The Non— Linear Case 

Let Q[i] be the field of Gauss rationals. For every positive integer number 
d G N, let Hd be the complex vector space of all homogeneous polynomials 
with complex coefficients of degree d. Namely, 

Hd ■= {/ G C[Xo, . . . , X„] : f homogeneous, deg(/) = d}. 

For every list of degrees (d) := (di, . . . , dn) € N", let 'H(rf) be the complex 
vector space given by the following identity : 

^(d) '■= X • • • X Hd^. 

The space "Hi^d) is the space of systems of complex homogeneous polynomials 
F := (/i, . . . , /„) such that deg(/i) = di. 

For every list of degrees (d), let D(^d) be the maximum D(^d) '■= maxjdi, . . . , d„} 
and let be the Bezout number of (d), namely V^^d) '■= TlILi ^i- 
Let Nd be the dimension of TCd- We have 

For every system F £ 'H^d)-, let ^norm{F) be the normalized condition num- 
ber introduced in |^9|, |3^] (cf. also |^ and Subsection below). This is a 
projective invariant (i.e. Hnorm{F) depends only on the class defined by F 
in the complex projective space F('H(rf))). 
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Let N be the (complex) dimension of the projective space F{Ti.(^ci))- Observe 
that the following equality holds : 



vj=i / \j=i 



For every list of degrees (d) := {di, . . . , dn), we denote by 'H(rf)(Q[i]) the set 
of systems of homogeneous polynomials F := such that fi G 

Q[i][Xo,...,X„]. 

Let P('H(rf)(Q[i])) be the projective space defined by the Q[i]— vector space 
'^(d)(QW)- In this case the (unitarily invariant) bit length of a projective 
point F £ P('H(^)(Q[i])) is defined as the logarithm of the unitarily invariant 
height of F (cf. Subsection |4.4| below). Then, the following statement holds : 



Theorem 3 (Main Theorem for the Non Linear Case) With the 
same notations and assumptions as above, there is a positive universal con- 
stant C2 G M, C2 < 20 such that the following holds : 

Let (d) be a list of degrees and let h,w G M be two positive real numbers. 
Assume that w > 1 and h is big enough to satisfy the following inequality : 

h > C2iV^ log N + log w. 

Then, for every system of polynomial equations F £ P('H(^)(Q[i])) of (unitar- 
ily invariant) bit length at most h, the condition number ^norm{F) satisfies 

linorm{F) < {w€.[{d)]f'\ (2) 

with probability at least 

2 

1--, 
w 

where 

n 

<t[{d)] := n^{n + 1)A^(A^ - 1) J] 

i=l 

The reader himself may rewrite this statement for the quadratic case (i.e. 

(d) = (2, 2, ... , 2)). Then he will realize what we mean. 

This statement is an immediate consequence of the following more technical 

(but more complete) statement. Li order to state this new Theorem, let us 

introduce some more notations : 

Let Sd be the set of multi-indices 

Sd := {(mo, •••,/««)£ N"+^ : |/i| = /io + . . . + Mn < 4- 
We introduce a well-ordering in the set of multi-indices Sd as a bijection 
(/) : Srf ^ {i G N : 1 < 1 < iVrf}. 
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Using this well-ordering on S^, we define the diagonal matrix A^^ in the 
following terms : 

Ad-.-- 

\ \ Lin . . . U'r, I I 

l<0(/iO,...,/in)<A'd 

where 

d \ / dl 




For every list of degrees (d) := {di, . . . , dn), let A(^) be the diagonal matrix 
given as the diagonal sum of A^-^ , ■ ■ ■ , Ad„ . Namely, 

A(rf) := Ad, e...eAd„. (3) 

Theorem 4 With the same notations as above, let H,e £M be two positive 
real numbers. Assume that H > 1 and e > 0. Let (d) := {di, . . . ,dn) be a 
list of degrees. 

Then, for every system F € P('H(rf)(Q[i])) of unitarily invariant height at 
most H the condition number ^norm{F) satisfies 

(F) > 1/e 

with probability at most 



C{2N + 2)~^H-g{N) 
where ^[{d)] is the constant introduced in Theorem [J above and 

q{2N+2) I 
^^^^ det(A(rf))C(2iV + l)i^2iV+2 ^ K2N+2 ^ ^' 

and 

Tn ■■= max{8,4(Z)(,) + i)}4[^^+3{n+2)2(iv+i)] ^ 
1.3 An application 

We apply the estimates of Theorem ^ above to estimate the probability 
distribution of the precision required to write down affine approximate zeros 
in Q[i]" of systems of multivariate polynomial equations of bounded input 
length. The reader may look at Section ^ for detailed definitions. For every 
z G the precision is defined as the bit length of a denominator of 

z. For every system of polynomial equations F £ 'P(d)(Q[^])) the precision 
Pr[F) of F is defined as the maximum precision required to write down 
approximate zeros in Q[i]" for every actual zero of F. Then, we show the 
following statement. 
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Corollary 5 With the same notations as above, there is a universal con- 
stant C3 > ( C3 < 20) such that the following holds : 

Let (d) := (di, . . . ,c?„) he a list of degrees and let h,w £ W be two positive 
real numbers. Assume that w > 1 and h is big enough to satisfy the following 
inequality : 

h > c^N^ {log N + log D^^^)) + log w, 

Then, for every system of polynomial equations F G P('P(^)(Q[i])) of bit 
length at most h, the required precision Pr{F) to write down an approximate 
zero of F satisfies 

Pr{F) < 0{n log2 D^^-^ + log^ w), (4) 
with probability at least 




1.4 Additional Comments 

As we have used different notions coming from different fields and different 
approaches, and we want to make our pages as readable as possible, we 
have included most notions and most basic facts in separate Sections. For 
instance, in Section § we introduce the basic notions and notations to deal 
with the Linear Algebra case, whereas in Section ^ we do the corresponding 
task for the Non-Linear case. Thus, we show Theorems |l] and || in Section ^ 
and Theorems |^ and § in Section ^. Finally, Section ^ is devoted to introduce 
the basic notions and notations to show Corollary ^ above. 
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2 Basic notions and notations 



2.1 Notations for the real projective space 

For every positive integer number m G N, let P^(M) be the real projec- 
tive space of dimension m. We denote by ttr : M"^+^ \ {0} — > Pm(lR) 
the canonical projection onto the projective space. Namely, given X := 
{xo,...,Xm) G \ {0} we denote by ttk{X) G Pm(M) the projective 

point whose homogeneous coordinates are given by the following identity : 

■kr{X) := {xq : xi : . . . : Xm) G P 



Let < •, • > : R™-"*"^ x — > M be the canonical Euclidean inner product 
in M™+i. Namely, given X := {xq, Xm),Y := (yo, • • • , ym) e M'^+i we 
define 

m 

<X,Y >:= ^Xipi. 

i=0 

For every X G M"*"^-*^ we denote its canonical Euclidean norm as ||X|| := 

(< x,x>)^/'^. 

This Euclidean inner product induces a natural Riemannian structure on 
Pto,(IR) which we denote by (Pm(I^), can). The Riemannian metric on 
(P„i(M), can) is denoted by (Ir. The metric (Ir is given by the following rule : 
given ttk{X) , ttk{Y) G Pm(IR), the Riemannian distance between Tr^{X) and 
t^rO^) is given by the following identity : 

<X,Y>\\ 



(iij(7rK(X),7rK(y)) := arccos- 



l^lllll^l 



We also introduce a Fubini-Study metric dps 

on Py^(M) from the previous 
Riemannian metric in the following terms : given 7r]^(X) , 7r]^(Y) G Pm(K), 
the Fubini-Study distance from ir^.iX) to Tr^{Y) is given by the following 
identity : 

dFs{'^R{X),7ru{Y)) := sin dR{TTR{X),-KR{Y)). 
For every V C P„(M), let V C M™+-'^ be the cone over V, namely 

F:=7r^i(F)U{0} CM'"+i. 
The Fubini-Study metric satisfies the following useful property. 

Lemma 6 Let V C Pto(IR) be a subset of the real projective space and let 

V := 7r^^{V) U {0} C 1^^"+^ be the corresponding cone over V . Then, for 
every point X G M"*"*"^ such that \\X\\ = 1, the following equality holds : 

d{X,V) := dFs{MX),V) = in f {dps MX), MY)) ■ MY) e V}, 

where d(X,V) := m/{||X — : Z G V} and \\ ■ \\ is the canonical Euclidean 
norm. 
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Let be the unit sphere in R"*+^, namely 

5™ := {X G R^+i : ||X|| = 1}. 

The canonical Euclidean inner product in also induces a Riemannian 

structure on 5™ which we denote by (S"",can). Let us also denote by 
PR '■ S"^ — > IPm(l^) the natural projection. Namely, 



For every measurable subset U C we denote by Vol{U) the standard 

Lebesgue measure of U, i.e. 



/ 

Jm. 



m+l 



Xifdxi ■ ■ ■ dx 



m+l) 



where xu '■ l^"*"*"^ — ^ M is the characteristic function associated to U. 
The canonical Riemannian structures on Pm(M) and S'^ respectively yield 
volume forms on every space. These volume forms lead to the following 
notions and notations. 

For every measurable subset U C S'^ we define the spherical volume of U 
in the following terms 



VolsiU) := [ xu{0)de, 



where d6 is the volume form in spherical coordinates. Observe that 

m+l 

VolsiS^ = 



r(^)' 

where F is the gamma function and xu • S"^ — ^ ^ is the characteristic 
function associated to U. 

Finally, for every measurable subset U C P„(R), its projective volume sat- 
isfies : 

VolpiU) := \volsipi\U)). 



2.2 Condition numbers for square matrices with real entries 

As observed in the Introduction, the usual condition numbers of Linear 

Algebra are naturally defined in the projective space P(A^„(M)) = P^2_;^(M). 

Here we recall these notions and some standard facts about them. 

Let S C A4„(R) be the algebraic variety of all n x n singular matrices. 

Namely, 

S ■.= {M e Mni^) : det(M) = 0}. 
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Let us identify A4„(M) ^ R" . The norm associated to the canonical Eu- 
chdean inner product in A^„(R) is usuahy called the Probenius (also Hilbert- 
Weil, also Schur) norm. Given a square n x n real matrix M £ A^„(M) we 
denote its Frobenius norm as HMHir. We also denote by dp{M,N) the 
Frobenius distance of two matrices M,N G A4„(R). We denote by ||M|| its 
norm as linear operator, namely 

||M|| := sup{\\Mx\\ : X G R" ||x|| = 1}. 

An outstanding statements concerning the norm ||M|| is the following state- 
ment due to C. Eckardt and G. Young. 

Theorem 7 For every M £ A^„(R), 

"^^ ^" " dpiM, 5)' 
where dp is the Frobenius distance. 

The usual condition number of a square matrix M G A^„(R) is given by the 
following identity : 

k{M) := ||M||||M-^||. 

A second (and useful) condition number /i was discussed in 
be defined as follows: for every M G A4„(R), we define 

/i(M) := ||M~i||||M||i.. 

As observed in ^, the following inequality holds for every matrix M G 

k{M) < ii{M), VM G 7W„(R). (5) 

For every matrix M G A^„(R), both condition numbers k{M) and f^{M) 
depend only on the projective point 7rR(M) G P(A1n(R)). Thus, we preserve 
the notation k and to denote the mappings /c,/i : P(A1„(R)) — > R. 
Namely, A;(7rR(M)) := k{M) and ;u(7rK(M)) := n{M). Let us denote by 
n(R)) given by the following identity : 

S:=MS\{0}). 

For every projective point 7rK(M) G P(A^„(R)) we define the distance to the 
algebraic variety of singular matrices S as 

piMM)) := dps [mm), ^) := infidpsMM), Z) : Z € S}. 

The following statement is an outstanding consequence of Eckardt & Young 
Theorem above. 



34, |. It can 
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Corollary 8 For every Z e P(A4n(M)), the following holds 

1 



k{Z) < fi{Z) 



P{Z) 



Proof- This fact is a consequence of Theorem ^ and Lemma |^ above. 
The ideas are essentially due to [^] (cf. also (P). We include a proof for 
didactical reasons only. 

The first inequality is Inequality |5| above. As for the equality, let M € 
Mn{^) be a square matrix such that 7rij(M) = Z. Let M' e 5"^"^ C 
Mn(M) be the square matrix given by : 

M' := „ \ M G MJM.). 

\\M\\f 

We obviously have 

k{M') = k{M) 

and ^m(M) =pr{M') = Z. 

Moreover, as ||M'||i? = 1, we conclude from Theorem ^ above that the 
following equality holds : 

fi{M) = fi{M') = \\M'-^\' - ^ 



dF{M',S) 



Now, observe that 5* is a cone. Namely, S = tt^^{S) U {0}. Then, from 
Lemma ^ above, we conclude : 

dF{M',S) = dps [MM'),s) = p{Z), 

and the statement follows. ■ 

For every e > 0, let 5(e) be the compact neighborhood of S in P(A^„(M)) 
given in the following terms : 

S{e) := {Z G nUnm ■■ dFs{Z,S) < e} = {Z e P(7W„(M)) : p{Z) < e}. 

The following statement is a mere consequence of Corollary ^ above. 

Corollary 9 

S{e) = 7T^HS{e)) U {0} = {M e : fi{M) > -J. 

The following Theorem can be found in pi]. 
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Theorem 10 With the same notations and assumptions as above, the fol- 
lowing holds for every e > : 

Volp{S{e)) < eVolp (P(7W„(M))) Tt''^. 

Sharper estimates can also be found in 0, |^] but they do not change es- 
sentially our forthcoming arguments. From this Theorem we may easily 
conclude the following statement. 

Let i?, e G M be two positive real numbers, e > 0. Let B{0,H) C A^„(M) 
be the closed ball of radius H centered at the origin with respect to the 
canonical Euclidean norm. We define the compact subset 

R{e, H) := B{0, H)f]§(^) = B{0, H) f| [7T^}{S{e)) U {0}) . (6) 

Corollary 11 Let H,e gM. be two positive real numbers, e > 0. With the 
previous notations and assumptions, the following inequality holds : 

Vol{R{e,H)) < en^/^Vols{S'''-^)H''\ 

Proof- Let us denote by R the compact set R(e,H). By integration in 
spherical coordinates we obtain the following equality. 

Vol{R):= [ XRdxi,i---dXn,n= [ [" XRirOy'^'-'drde. 
JR"^ JS"'^-^ Jo 

Now, let us observe that XR is homogeneous of degree (for < r < H). 
Namely, given r, r' G R, < r,r' < H, XRi^G) = XR{f'(^)y for ah 0. 
Then, since vr^(i?) = S{e) we conclude the following identity : 

Vol{R) := {^j^^^_XR{e)de^ ir^^^''''^') ^ ^Vols{pi\S{e))), 
and the statement follows from Theorem |l^ above. ■ 

2.3 Davenport's discrepancy Bounds 

Here we recall a beautiful statement due to H. Davenport (cf. |5|). Let 
m G N be a positive integer number. 

Given a measurable subset R C M™- of finite volume Vol{R) < oo, we denote 
by N(R) the number of points in R with integer coordinates. Namely, 

N{R) :=tt(i2n^™) • 
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It is well-known (cf. j2[l| ], for instance) that the following equality holds : 

where /i € M is a positive real number and hR := {hx : x € R}. Upper 
bounds for \\N{hR) — Vol{hR)\\ are usually called discrepancy bounds (cf. 



[10 1, |11|, and references therein for other approaches). Here we recall 
a method due to H. Davenport to obtain discrepancy bounds which by the 
way is very accurate for the problems we deal with. 

For every positive integer number £ G N, 1 < i < m, we denote by 2^"^^ the 
class of all subsets I of {i : 1 < i < m} such that lj(/) = i. 
For every subset / G xl™"^ , let us denote by the projection 

: M™ — > R^, 

given by 

^i{xi, . . .,Xm) ■■= {xi : i £ I). 

Let i? be a compact subset of M™". We denote by V{R,£) the following 
quantity : 

V{R,e):= ^ Vok{^j{R)), 

/ex*™' 

where Vole{^i{R)) denotes the Lebesgue of $/(i?) as subset of M^. We 
define V{R, 0) = 1 by convention. 

For every subset U C R"^, we denote by (5q{U) the number of connected 
components of U (i.e. the 0— th Betti number). 

Let / € X^™"^ be given, and let J <^ I, such that {](/ \ J) = 1. For every 
point A := (ctj : i £ J) £ M^""^, let rj{A) be the real line in given by 
the following identity : 

rj{A) := {{xi -.i £ I) : Xk = ak,yk £ J}. 

Let R C be a compact set. We define h{R,£) as the following number : 

h{R,i) := sup{(3o{^i{R)nrjiA)) : I £4"'\j C I,UI\J) = l,A£R'-^}, 

where Pq denotes the 0-th Betti number (i.e. the number of connected 
components. 

We finally define h{R) := sup{h{R,£) : 1 < ^ < m - 1}. 
These notations stated, the following statement holds. 

Theorem 12 (Davenport's discrepancy bound) ^ Let R C M™ a com- 
pact set. Let N{R) be the number of points in R with integer coordinates. 
Then, the following holds : 

m—l 

\N{R) - Vol{R)\ < h{Rr-^V{R,i). 

1=0 
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2.4 Northcott— Schmidt's height, bit length and visible points 
from the origin 

Let A C M'^+i be a lattice. A non-zero point X E A \ {0} is said to be 
visible from the origin (or simply visible) if there is no point of A in the 
segment [0, X] C M"*+^ between the origin and X. 

For every non-zero real number p G M\{0}, let Z'""^-^(p) be the lattice given 
by the following identity : 

Z"^+i(p) :={pX : XgZ"*+1}. 

It is easy to prove that the set of visible points in (/?) is the set of all 
those points pX := {pxo, . . . ,pxm) such that the greatest common divisor 
gcd{xo, . . .,Xm) = 1. 

For instance, let A^„(Z) be the set of all n x n matrices with integer entries. 
The set A^„(Z) is a lattice in A4n(IR) — IR" • For every non-zero real 
number p G M, let A^„(Z)(/9) be the lattice in A^„(M) given by the following 
identity : 

X„(Z)(p) := {pM : Me 

The visible points of A^„(Z)(p) are exactly those matrices 
pM := {pxi,j)i<i,j<n such that the greatest common divisor 
gcd{xij : I <i,j <n) = 1. 

Visible points in a lattice Z'^'^^ (p) and projective points are closely related. 
In order to illustrate this relation, let us introduce the Schmidt height of a 
projective point. Let P C N be the class of all prime numbers. For every 
p e P, let II "lip : Q — > M+ be the non-archimedean p-adic absolute value. 

Definition 13 Let TrM.{X) := (xq : ... : Xm) G ^miQ) o, projective point 
whose homogeneous coordinates are rationals, i.e. assume that 
X = {xo,...,Xm) G Q"""*"^ \ {0}. We define the Northcott-Schmidt abso- 
lute height of ttm.{X) as the following infinite product : 

H{Trm{X)) := \\X\\ JJ max{||xi||p : < i < m}. 
peP 

This notion is well-defined because of WeiFs product formula. Let us con- 
sider now the projective space P(A^„(Q)) define by the vector space of all 
n X n matrices with rational entries AiniQ)- Given a non— zero matrix 
M := (xij) e A^„(Q), the Northcott-Schmidt absolute height of 7r]R(M) is 
given by the following identity : 

i7(7rM(M)) := ||M||f JJ "^^iW^idWp 1 < < 
peP 

where ||M||i^ is the Frobenius norm introduced above. 
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Let pX := (pxQ, . . . , pxm) £ be a visible point in the lattice 

Z"^~^^{p). Let 7rR(pX) = tt]r{X) G IPm(Q) the corresponding projective 
point. Then, it is easy to show the following Lemma. 

Lemma 14 With the same notations and assumptions as above, the follow- 
ing properties hold : 

i) Let pX := (pxo, • ■ • ,pxm) £ be a point in the lattice 7/^'^^{p). 

Then, pX is a visible point from the origin if and only if 

\\X\\ = min{\\Y\\ : pY G Z"'+\p),Mpy) = M^)}- 

a) Let pX € Z'""*"^(p) be a visible point. Then, 

\\X\\ = H{ttm{pX)). 

Hi) Let V C P(A4„(M)) be a subset of a real projective space and let V := 
7r~^{V) U {0} C Mn{^) the cone over V. Let B{0,H) be the closed 
ball in R"*+^ of radius H centered at the origin. Then, the following 
holds : 

• For every i^-m.{X) G F PI Pm(Q) there is a visible point Y e V D 
Z'^+i such that irmiY) = TrR{X) and \\Y\\ = H{TrM{X)). 

• The number of visible points of 'L^^^ that lie in V Ci B(0,H) 
is two times the number of projective points in V Ci Prn(Q) of 
Northcott- Schmidt absolute height at m,ost H. 

In fact, given a projective point 7rR{X) G Pm(Q) and given a visible point 
Y e Z"'+^ such that m{Y) = 7ru{X), the two visible points in Z"*+^ associ- 
ated to 7Tr{X) are {Y, -Y} C Z"*+i. 

Given a projective point Tr^{X) G Pm(Q); we may represent this point by 
some visible point Y G Z™+^ such that 7rE(y) = TT^iX). Thus the bit length 
of irM.{X) is the number of tapes cells required to encode y in a Turing 
machine. As Y G Z"*+^, this quantity is essentially equal to the logarithm 
of the norm of y, i.e. log ||y||. This justifies the following Definition. 

Definition 15 Given a projective point '7rM.{X) G Pm(Q); we define its bit 
length as the logarithm of its Northcott- Schmidt absolute height. Namely, 

bl{MX)) = log2H{Tr^{X)). 
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3 The Linear Algebra case 

3.1 Davenport's discrepancy in the Linear Algebra case 

In this Subsection we apply Davenport's strategy (as introduced in Subsec- 



tion 2.3 above) to count the number of integer points in the region R{e, H) 



introduced in Equation (|6|) of Subsection 2^ above. Namely, given two pos- 
itive real numbers H,e G M, let R{e,H) C A^„(M) be the compact subset 
given by the following identity : 

R{e, H) := B{0, H)f]§(^) = B{0, H) f| [7T^}{S{e)) U {0} 

Proposition 16 With the previous notations and assumptions, let H > \ 
he a real number and let iV(e, H) he the number of matrices with integer 
entries in R{e,H), namely 

N{e,H) ■.= ^(^R{e,H)f~]Mnm 

Then, the following inequality holds : 

\N{e, H) - Vol {R{e, H)) \ < r„e("')i7"'-^ 

where ©("^) is the constants introduced in Equation |I| of the Introduction 
and 

T„:=(2max{4,n})(2"'+^"'), 
is the constant introduced in Theorem |J. 

In particular, from Corollary |71| we conclude the following inequality 
N{e,H) < eni/Vo/s(5"'"^)i^"' + i;6^"'^i^"'"^ 

Proof- First of all, observe that R{£, H) is the set of all matrices M G 
AinO^) that satisfy the following properties : 

• ||M|||, < 

• dFsMM),S) <£. 

Now, let M' G A^„(]R) be the square matrix given by the following identity : 

M' := jT-^M. 

\\M\\f 

We obviously have, = 1 and ttk{M') = 7r/j(M). From Proposition ^ 

above, we conclude : 

dpsMM), S) = dpsiMM'), S) = dF{M', S). 

Then, the previous two properties that define R{£, H) may be replaced by 
the next two ones : For every M G Mn(M), M G R{e,H) if and only the 
following two properties hold : 
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\\M\\l < H'^, 

3N £ S, such that \\M - \\M\\fN\\f < eWMllp- 



Then, with the notations of Subsection 2.3 above, for every integer number i, 

1 < £ < n?, let Tf' be the class of all subsets I of : 1 < J < n} such 

that Jj(/) = i. For every I G T^^"'\ let <!>i : A4„(M) — > the projection 
defined by I in the following terms : 



{{xi,j)i<i,j<n) ■■= {xr,s ■■ {r,s) G /) G M^. 
Observe that for every / G 2^ the following holds : 

where B£{0,H) is the closed ball of radius H centered at the origin of M^. 
Then ifVol{^i{R{e, H)) denotes the ^-volume of ^i{R{e, H)), the following 
inequality holds : 

Vol{^i{R{£,H))) < Vol{Bi{0,H)) = H^Ke. 



With the same notations as in Subsection |2.3| and Theorem 12 above, the 
following inequality holds : 

V{R{e,H),£)= J2 Vol{^jiR{e,H)))<H'(^'^^Ke. 



The second invariant in Davenport's Theorem 12 is given in the following 
terms. Let £ G N be a positive integer number such that 1 < £ < n?. Let 

/ G I^"' ■* be a subset of : 1 < i,i < JT-} of cardinal £. Let J C / a 

subset such that ^{I \ J) = 1. Let A = {ci{i,j)){i,j)eJ ^ ^^^^ be any affine 
point and let rj{A) be the real line in given by 

rj{A) := {{xij : G I) G : Xij = Oij, G J}. 

We define h{R{£,H),£) as the maximum number of connected components 
of the sets given by 

^i{R{e,H))nrj{A), 

where I C T^^'^ J C / such that tt(/ \ J) = 1 and yl G M^-^. 
With the same notations, let us introduce the real algebraic subset 
V{e, H, I, J, A) C M?^ given by the following property : a point 
{M,N,ti,t2,t3,U) G MnO^f X ^ m2"'+4 belongs to V{e,H,I,J,A) if 
and only if the following properties hold : 

i) Xij = Uij, for all (z, j) G J, 
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ii) H'-\\M\\l = tl 

iii) det(7V) = 0, 

iv) IIMIII. = 



V) t2-4 = 0, 

vi) \\M\\le^-\\M-t2N\\l = tl 

We observe that the maximum of the degrees involved in these equations is 
max{4, n}. In other words, the polynomials of maximal degree occur either 
in the equation iii) or in equation vi). Now we may apply the Milnor- 
Thom-Oleinik-Petrovsky upper bounds for the sum of the Betti numbers of 



real algebraic varieties (of. ||2^, ||2J]). These upper bounds yield 

the following inequality : 

A (Vie, H, I, J,A))< (2 max{4, n})^"'^^ . (7) 

On the other hand, let ^ : Mn(M)^ x — > be the projection 

given hy ^{M,N,ti,t2,t3,t4) = M. Then, the following identity holds : 

CD, (M/ {V{e, H, I, J, A))) = ^i{R{e, H)) n rj{A). (8) 

In order to see this identity, let us observe that Equation ii) above means 
11-^11 F ^ H, Equation iii) means that A'^ is a singular matrix and. Equations 
iv),v) and vi) mean that \\M — \\M\\pN\\ < \\M\\pe. 

Finally, as ^ and are continuous mappings, we conclude from Identity 
(^) and Inequality (|^ above the following inequality : 

/3o{^i{R{e,H))nrj{A)) < l3o{V{e,H,I,J,A)) < (2max{4,n})2"'+^ 

In particular, for every i ^ 'N, 1 < £ < n? , the following inequality holds : 

h{R{e,H),i) < (2max{3,n})2"'+^ 

Finally, as if > 1 we conclude from Theorem |l2| the following inequality : 

\N{R{e,H)) - Vol{R{e,H))\ < H""'-^ (2max{3,n})(2"*+^"') 



Proposition 17 With the same notations and assumptions as in Propo- 
sition above, let H > 1 be a positive real number. Let N{1,H) be the 
number of matrices with integer entries in R{1,H). Namely, 

N{1,H) :=s(s(0,i/)f|7W„(Z)) . 
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Then, the following holds : 

\N{l,H) -Vol{B{0,H))\ < i/"'~^6("'). 
In particular, the following inequality also holds : 



Vols{S' 



Proof- The proof follows the same arguments as in Proposition ^ above, 
noting that h{B{0, H),i) is always 1. ■ 



3.2 On ill— conditioned matrices that are visible from the ori- 
gin. 

Proposition 18 Let H,e £ M be two positive real numbers, H > 1. Let 
M{e, H) be the number of points mP(7W„(Q)) of absolute Northcott- Schmidt 
height at most H which belong to S{e). Then, the following holds 



where 

r„6("') , Vol{R{e,l)) 
2C(n2 - 1) 



L{e,n) ■= \ ^^^^ — + 



and Tn and K^^ ^ are as in Proposition [7^ above. 

In particular, the following inequality is a consequence of Corollary 11 : 



where 



Proof.- This proof is widely inspired by the Proof of Theorem 459 of |14]. 
Assume that e > is fixed throughout this Proof. Let V{e) C 7W„(M) = M"^ 
be the compact subset given by the following identity : 

V{e) :=S(0,l)n(^-'(S(£))U{0}). 

From Corollary we easily conclude that for every positive real number 
t G M+ 

tV{e) := {tX : X £ V{e)} = R{e,t), 
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where R{e,t) is the compact set introduced in Subsection 3.1 above. In 
particular, for every positive real number t £ M+ the following identity 
holds : 

r\ol{V{e)) ■.= Vol{R{s,t)) (9) 

For every positive real number p £ M+, let us denote by Ap the finite subset 
of i?(0, 1) given by the following equality 



Ap:=Vie)f]Mnmp)\{0}, 



where is the lattice in 7W„(M) introduced in Subsection 2.4 . 

Namely, 

Mnmip) := {pM : Me 

Let g{p) be the number of points in Ap. Observe that the following equality 
holds : 

g{p) := tt(^p) = tt {R{e, p'') f] A1„(Z)) - 1 = N{e, p-') - 1. (10) 

Let f{p) be the number of points in Ap which are visible from the origin. 
As in the Proof of Theorem 459 of |14] we may easily conclude that 

oo 

9{p) ■=^f{mp). 
1=1 



From Mobius inversion formula (see Theorem 270 in |14|, for instance) it 
follows that 

oo 

fip) ■= P('^)9{'mp), 

m=l 

where p is Mobius function. Observe that Riemann's zeta function satisfies 
the following identity for every s > 1 (cf . Theorem 287 of |0] , for instance) : 



^-E^- (11) 



C{s) ' ^ m 
So 

""'^f"' - ^tS^ = £ {Cnpr'simp) - Vom.))) ^ (12) 

^ ' 771=1 

Let Np be the set of integer numbers given by the following identity 

Np := {m G N : 1 < m < p'^}. 

Let us denote by Sp the sum of the terms of the left hand series of Equation 
(30) whose indices are in Np. Namely, 



^p--= E ^({mpfgM-VoliVie))). 

meNp 
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Now, assume that m E Np (i.e. mp < 1). Then, replacing (mp) ^ by H, we 
conclude from Identities (^) and (|l^) above the following identity : 

\{mpfg{mp) - Vol{V{e))\ = -^\N{e,H) - (1 + Vol{R{e, H)))\ (13) 



From Proposition 16 above, we have : 



1 T ) 

\N{e,H)-Vol{Rie,H))\< " 



Replacing back H by (mp)"^, Equation (|l3| ) becomes the following inequal- 
ity : 

|(mp)"'g(mp) - yo;(y(e))| < r„©("')mp + (mp)"' (14) 



We conclude : 



E (;$^''^»e("') + M™v"'). (15) 



meNp 

Then, 

S.<^^^ + ^E^irn)p j<^^^^p . (16) 

On the other hand, let be the set of positive integer numbers which are 
not in Np. Namely, 

:= {m G N : mp > 1}. 
Now, assume m G and X G Amp- Then, we have X = (■mp)Y, where 

2 

y G . In particular, the following inequality holds : 

||y|| = (mp)||y|| > mp > 1. 

In other words, for every m G N'^, g{mp) = 0. Hence, let us define S'p as the 
sum of the terms in the left hand series of Equation (12) above. Namely : 



S'p--= E ^\{mpf9{rnp)-Vol{V{e))\. 

meN' 



We have : 

S'p<Yl ^\Vol(V{e))\<V{e) I ^ | <pVol{V^)) (17) 

meN^ \m>p-l / 



Now we combine Inequalities ( |16D and (17) to conclude : 

^^^^ C(n2) I -'^'^ + '^/-- C(n2-1) +^ +pFo/(y(e)). 
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Finally, replacing p ^ by in this Equation and using the identity described 
in Equation (P) we conclude : 

) ^) — '-cR^^ +^ + H • 

Noting that f{H^^) = l/2J\f{e, H) (cf. Lemma 14 above) the Proposition 
is achieved. ■ 

Proposition 19 Let H gM. be a positive real number, H > 1. Let H) 
be the number of points in P(A4„(Q)) of Northcott-Schmidt absolute height 
at most H . Then, the following holds 



whe 



re 



■= 2C(n^ - 1) + ^ 



2 



where Tn and Kf' are as in Proposition 16 above, and ( is Riemann's zeta 
function. In particular, the following inequality also holds : 



Proof- This is step by step the same proof as that of Proposition 18 using 
Proposition ^ instead of Proposition 16. ■ 



3.3 Proof of Theorem |l] 

Proof of Theorem ||.- First of all, observe that the probability that a random 
choice of a system of linear equations M G P(7W„(Q)) of Nosthcott-Schmidt 
height at most H satisfies ii{M) > 1/e is given by the following quotient : 

M{e,H) 



Let A be the quantity : 

^ K„2 



2C(n2)' 

From Propositions O and |l^ we conclude : 



Af{e,H) ^ en'>I^AH + L'{e,n) + ^^ 



N{l,H)- 
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Hence, we have 



< 



Finally, as > 1 we conclude : 



L{l,n) _^ 1 



N{e,H) 



N{l,H) 



en 



5/2 



L{l,n) _^ 1 



if 



where 



and 



B{e,n) :=en^/2 

C(n): 



L(l,n) ^ 1 



A 2A 



L'(e,n) 1 
+ . + 



^ 2A 



4 Condition Numbers for Systems of Multivariate 
Polynomial Equations 

4.1 Two Hermitian inner products in 7i((i). 

For the sake of readability, we recall here some of the notations stated in 
Subsection of the Introduction. We define Hd as the set of all complex 
homogeneous polynomials in n + 1 variables of degree d. Namely, 

Hd ■= {/ G C[Xo, ...,Xn] : / homogeneous, deg(/) = d}. 

We define the standard Hermitian inner product in by identifying Hd = 
C^'^ , where A^^ is the number of coefficients of a generic homogeneous poly- 
nomial / G C[Xo, . . . , Xn] of degree d. We have 



Nd-.-- 



d + n 
n 



The standard Hermitian inner product < • , • > : Hd x Hd — > C is given 
in the following terms : Given f,gG Hd we define < f, g > as 

< f,g >■= J2 "a^^M' 

where fi := {fiQ, . . . , fin) G N"^"*^ is a multi-index, \fi\ := /xq + • • • + /^n, 
f:=Y^ a.Xi;^ • • • X^, 5 := E «M^o" " " " ^n^" 

\fj.\=d \fj.\=d 
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and stands for the complex conjugate of b^. 

We introduce a well-ordering in the set of multi-indices Sd ■= {{fJ-o, • • • , /^n) '■ 
Ho + ■■■ + fin = d} as a bijection 

^ :Sd^{l<i<Nd}. 

The reader may assume that ip is given by the lexicographic order on S^- 
Using this ordering we define the diagonal matrix in the following terms : 




where 



^<'P{lJ-0,---,IJ-n)<Na 

d \ / dl 



As in Q, we also define an Hermitian inner product < ■,■ >d ■ Hd x Hd — > 
C in the following terms. 

< f,9 >d-=< ^df,^dg > ■ 

For every list (d) := {di,. . . , dn) S N" of degrees let Hi^d) = -f^di x • • • x Hd^ 
be the space of sequences F := (/i,---,/™) of homogeneous polynomials 
such that fi € if^. , 1 < i < n. We may also see 'H(^d) as the space of all 
polynomial mappings F := (/i, . . . , /„) : C"^^ — > C", such that fi G Hd^, 
for all i. 

We may extend the previous Hermitian inner products to the product space 
'H(^d) ill the following terms. Given F := (/i,...,/n) £ "^(d) and G := 
{gi, . . . , Qn) € T~(-{d)i define the canonical Hermitian inner product on 
Ti-(^d) as < •, • > : Ti-i^d) ^ "^(d) — * ^ in the following terms : 



< F,G >:=^< fi,gi > 



i=l 



For every F G 'H(d), we denote by ||-F|| the norm defined by this canonical 
Hermitian inner product, i.e. ||-F|| := (< F,F ■ 

On the other hand, let A(^d) be the diagonal matrix given as the diagonal 
sum of , . . . , Ad„ . Namely, 

A(rf) := Ad, (B ■ ■ ■ (B Ad„. 

Then we also introduce the Hermitian inner product used in ||3| in the fol- 
lowing terms. Given F,G £ 'H(^d)j 'we define < F,G >aS K as 

<F,G >a: = < A(rf)F, A^d)G >G M. 
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For every F E 'H(d), we denote by ||-F||a the norm defined by this canonical 
Hermitian inner product, i.e. ||-F||a := (< F-iF >a)2. 

Let U{n + 1) be the unitary group of isometrics of C""*"^. We may define the 
action of the group U{n + 1) on 'H^d) iii the following terms. Given F G 'H(^) 
and given cr G C/(n + 1), we define o"(F) € 'H(^) as the unique polynomial 
mapping in 'H^d) that satisfies the following identity. 

(j{F){x) := F{a~\x)), Vx G C"+^ 

Then, the following statements holds 

Theorem 20 ^ This Hermitian inner product < •, • >a on Hi^d) ^-^ uni- 
tarily invariant. In other words, 

< a{F),a{G) >a=< F,G >a, 
for every a £ U{n + 1) and for every F,G £ T~(-(d) ■ 



4.2 Two Riemannian structures on F(H(^d))- 

Let F{Tl(^d)) be the complex projective space defined by the complex vector 
space Tii^d)- Hence, S N is the complex dimension of F{TC(^d))- Recall that 
the real dimension of ¥{7i(^d)) is 2A^ and that N is given by the following 
identity : 

The well-order on the monomial exponents introduced in in Subsection 
of the Introduction, yields an identification between 'H(^d) C^^^. Ac- 
cordingly, we may also identify P('H(^)) and PAr(C). 

We may introduce two Riemannian metrics in PAr(C) according to the Her- 
mitian inner products on C^"*"^ discussed in Subsection 4.1 above. These 



two Riemannian metrics yield two Riemannian structures in Fjv(C). We 
denote by (PAr(C),con) the Riemannian structure defined by the canonical 
Hermitian inner product < •,• > : 'H(^d) ^ ^(d) — *■ C. We denote by 
{¥ ]\f (C) , uni) the Riemannian structure defined by the (unitarily invariant) 
Hermitian inner product < •,• >a : T~(-(d) ^ ^{d) — * C introduced in 



Subsection |4J| above. 

Finally, let us denote by vrc : C^+^ \ {(0, . . . , 0)} — > PAr(C) the canonical 
projection given by 

7rc(2;o, ... ,zn) := {zo : zi : . .. : zn) £ lP(W(d)), 

where {zq : zi : . . . : xn) are the homogeneous coordinates of the projective 
point defined by (zq, . . . , zn) S C^"^^. 
Then the following Lemma holds : 
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Lemma 21 The following mapping is an isometry between (P('H((^)), can) 
and i^{'H(^fj^)^uni) : 



: (PAr(C),can) — > {¥]s{C),uni), 

where 

A"^(zo : zi : . . . : zn) := 7rc(A^^^(zo, ^^i, • • -.zn))), 
and AT^^j G GL{N is the inverse of the regular matrix A(^) introduced 



in Subsection 4-1 above. 



We denote by A the inverse of the isometry A ^ introduced above. Observe 
that for every (zq : zi : . . . : zjy) G Fiy{C) the fohowing holds : 

A(zo : zi : . . . : zn) := 7rc{/^{d){zo, zi, . . .,zn))). 

We may associate two volume forms to Pjv(C) accordingly to the corre- 
sponding Riemannian structure. For every subset B C P7v(C) we denote 
as VolcaniB) the volume of B with respect to the Riemannian structure 
{Fj^{C), can) and we denote as VoluniiB) the corresponding volume of B 
with respect to the Riemannian structure {F ]\j (C) , uni) . From Lemma ^ 
above the following holds for every subset B C Fj\f{C) : 

Volcan{A{B)) = VoluniiB). (18) 

Now let S"^^"*"^ C ]R^^+^ be the unit sphere. The unit sphere has also a 
canonical Riemannian structure defined by the canonical inner product in 
-[^2Af+2_ Identifying R^^"*"^ = C^"*"^ we may also consider the canonical 
projection 

p := TTc I52JV+1 : (5^^"*"^, can) — > (P7v(C), can). 
The following Proposition is a well-known fact in Riemannian Geometry. 



Proposition 22 The mapping p is a Riemannian submersion. In fact, the 
following holds : 

This Proposition yields the following identity for every subset B C P;y(C) : 

Vols{j>-^{B)) = 27rVokan{B), (19) 

where Vols{p^^{B)) denotes the volume of p^^(5) with respect to the canon- 
ical volume form (spherical volume) in 5^^"*"^. 

Together with the Riemannian structures we have the corresponding Fubiny- 
Study metrics in P7v(C). We reproduce here this metric for the Riemannian 
structure {Fiy{C),um). It is similarly defined in the other case. Recall 
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that the Riemannian metric in (Pjv(C),uni) induces the distance function 
: PAr(C) X PAr(C) — > R+ given by the following identity : 



dR{7rc{X),7rc{Y)) := arccos 



Re<X,Y >A 



ll^llAniA 

where Re < X,Y >a is the "real part" of the complex number < X,Y >a£ 
C. 

We define the Fubini-Study metric associated to the Riemannian structure 
{FN{C),uni) as the mapping dps ■ ^n{C) x PAf(C) — > R+ given by the 
following identity : 

dFs{x,y) := sin dR{x,y), Va;,y G PAf(C). 

Observe that the Fubini-Study metric satisfies the following property. 

Proposition 23 Let B C Pjv(C^) be a subset and B := ir;^^ (B) U {0} . Let 
X G C^^^ be a point such that 

||X||a = 1. 
Then, the following identity holds : 

inf{\\X-YU ■■ Y &TT^\B)VJ {{)}} = dFs{MX),B) 
where dFs{'^c{X),B) := inf {dFs{'^c{X),y) : y ^ B}. 

An analogous statement holds for the corresponding canonical Riemannian 
structure on 



4.3 Solution and Discriminant Varieties 

We may define the solution variety 

V := {(F, C) G P(W(d)) X P(C"+i) : F(C) = G C"}. 

Proposition 24 ^ With the previous notations and assumptions, V is a 
smooth connected projective subvariety o/P('H(rf)) x P(C"^^) of codimension 
n. 

Let us define the projection tti : V — > P('H(rf)) and for every (FX) G V, 
let 

Z)7ri(F,C) -.T^f^qV ^TFF{n(^a)), 

the corresponding tangent mapping at {F,(). We say that (FX) G 1^ is 
ill-conditioned if Dtti(F,() is not an embedding. We define the set of ill- 
conditioned inputs T,' <ZV hy the following identity 

J:' :={{F,C)eV : det Z)7ri(F, C) = 0}. 
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Then, we define the discriminant variety as the algebraic variety 



S := 7ri(S') := 7ri({(F,C) G V : det Z)7ri(F, C) = 0}) C P(W(rf)). 

Let TT2 '■ V — > P(C"^^) be the projection onto the second factor. Namely 
7r2(F,C):=Cforall(F,C)GK 

Given ( G P(C""*"^), we define := tt2^{C)- Namely, is the algebraic set 
of all systems of polynomial equations F that vanishes at i.e. 

:= {{F,x) e P(?^(rf)) X P(C"+i) ■.{F,x)GVAx = (}■ 

Given (F, G V, we define the fiber distance to the discriminant variety S 
in the following terms : 

p{FX):=inf{dFs{F,G) : (G,C)GS'f|y}. 

Lemma 25 /Q/ The function p : V — > M is invariant under the action 
of U{n + 1) on V. Namely, given {F,Q G V and given a G U{n + 1), the 
following equality holds : 

p{a{F),a{C)) = p{F,C). 

We also define the fiber distance of system F G V C P('H(^)) to the discrim- 
inant variety T, as 

p{F) :=min{p{F,C) : i^(C) = 0}. 
Obviously, p{F) = when F G S. 

On the other hand. Let F := (fi, . . . , fn) G 'H(rf) be a system of homogeneous 
polynomials. Let F : C""*"^ — > C" be the regular mapping between affine 
spaces defined by F. For every z G C"+^ let DF{z) : C"+^ — > be the 
tangent mapping defined by the jacobian of F at z. For every G C""*"^, let 
:= {li; g C"+^ : < w,C >= 0} the orthogonal complement of ( G C"+^ 
with respect to the canonical Hermitian inner product. 
Let us denote by Diag{\\C\\'^'~^dy'^) the square diagonal matrix given by the 
following identity : 

[wcr^M^' ■■■ 

D^ag{\\Cr'-'dy') := : •.. : 

V ••• iicr"-^dy', 

Let F G 7Y(rf) be a system of polynomial equations and let C G Pn(C) be a 
non-singular zero (i.e. -F(C) = and the linear mapping DF{(^) : — > C" 
is a non-singular linear mapping). The normalized condition number of F 
at C is define by the following identity : 

/x™(F,C) := \\F\\a\\DF{C) \t' Diag{\\Ct''dy^)\\, 
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where ||-F||a denotes the norm of F with respect to the unitarily invariant 
Hermitian inner product in 'H(rf) introduced in Subsection iA_ above. We 
define the normahzed condition number fJ-normiF') the maximum : 

fLnorm{F) := max{//™(F, C) : C e Pn(C),F(C) = 0}. 

In these pages, we are not going to discuss the relevance of this notion 
both for the behaviour of the projective Newton's method and for the com- 
plexity of homotopy continuation methods to solve systems of multivariate 
polynomial equations. The reader interested may follow the comprehensive 
discussion in Chapter 14 of |^] and the references therein. The following cen- 
tral statement explains one of the main features of this normalized condition 
number. 

Theorem 26 (Condition Number Theorem, [^) With the previous no- 
tations and assumptions, for every system F G ¥{7i^^^), the following holds : 

1 



Pin 

In fact, the Condition Number Theorem above gives a more precise state- 
ment that we are going to use in the sequel. This is explicitly given in the 
following statement. 

Proposition 27 ^ Let F := (/i, . . . , /„) € 'H(rf) be a system of homoge- 
neous polynomials and C € Pn(C). Assume C) £ ^ '^''^d ||-F||a = 1- Let 
a & U{n-\- 1) be a unitary transformation such that ct((") := (1,0, ... ,0) = 
eo S C""*"^ . Then, the following identity holds : 

p{F,C) := dF{Da{F){eo)Diag{d-'^^),S), 

where 

i) S A^„(C) is the algebraic variety of singular matrices, 
a) Da{F){eo) G A4„(C) is the jacobian matrix : 

0.(F)M:=(^(e„)) , 

— 1/2 

in) Diag{d- ) is the diagonal matrix given by the following identity : 

(d-"' ... 



Diag{d;'/') :-- 



\ 
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iv) dp is the Frobenius distance introduced in Subsection 2.^ above, namely 
given X := (xjj) and Y := {yij) G A^n(C), we define 

dF{X,Y):=\\X-Y\\p:= { ^ \xij - y^,J\^ 

y<i,j< 

For every e > we define the following tubular neighborhood of the dis- 
criminant variety S. 

:= {FGycP(7^(^)):p(F)<e}. 
Now, the following statement holds. 

Theorem 28 With the previous notations, the following holds : 

Voluniine)) 



Voluniin^id)) 



<eM(d)], 



where '^[{d)] := n^(n+ l)A^(A^ — l)P(rf) is the constant introduced in Theorem 
1^ above. 

For every positive real number H > 0, let Bj\{0,H) C C^^^ be the closed 
ball of radius H centered at the origin G C^^^ with respect to the unitarily 
invariant norm || • \\^. Namely, 

B^{0,H) ■.= {Z £C^+^ : \\Z\\a<H}. 

On the other hand, let B{0, H) C C^^^ the closed ball of radius H centered 
at the origin G C^^^ with respect to the canonical norm || • ||. The following 
equality holds : 

i?A(0,i/) = A-^5i?(0,i/), (20) 

where A^^^) : C^^"*^ — > C^"*""*^ is the complex vector space automorphism 
given by the regular matrix A(^). 

With the same notations as above, for every e > and for every positive 
real number H > 0, let R/\{e,H) be the compact set given by the following 
identity 

RA{e,H) := (7T^\j:{e))[j{0}) f~]BA{0, H). (21) 

The following statement is a consequence of Theorem ^ above : 

Corollary 29 With the previous notations and assumptions, the Lebesgue 
measure of Ra{£, H) satisfies the following inequality : 

Vol(RA{e,H)) < e'^H^^+^^l 



where 

_ 2vryo/...(F(H(,)))g[(d)] _ 
^^ ^^ (2iV + 2)det(A(,)) - det(A(,))''^^''^J- 



32 



Proof- For every H, let us define the compact subset V{e, H) given by the 
following identity : 

V{e,H) := [n^\A{J:ie)))[j{0})f]BiO,H), 

where 

A(S(e))):={A(x) : x G 
By integration in spherical coordinates we obtain 

TT2N+2 

Vol{Vie,H)) = ——Volsip-H^im))- 
Now by Identity (|l|) ab ove we conclude 

Vol{V{e,H)) = ____Fo/,,„(A(S(£))). 



Finally, from Identity (|18| ) above, we conclude : 

2Tr.ff2N+2 

Vol{V{e,H)) = ——-Vol^r^^im)■ 

On the other hand, the following identity holds : 

vr-i(A(S(£))) = A(,)(7r^H5](e)))- 
As Ba{0,H) = A^^^B{0,H), we conclude : 

V{e,H)=Af^^)RA{0,H), 

Thus, we conclude 

and the Proposition follows from Theorem above. 



4.4 Gauss integers, unitarily invariant height and C— visible 
points 

The theory of visible points (as used in previous pages) is not directly appli- 
cable to the study of points in a complex projective space. In this Subsection 
we shall see how to modify that theory of visible points to count points in a 
projective space IPm(QH), where Q[i] is the field of Gauss rationals. 
Let = {a + bi : a,b £ 7j} he the ring of Gauss integers and let 
be the group of units. Namely := {1,-1, where = —1. As 
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Z[i] is a factorial domain, we may also introduce the corresponding theory 
of visible points in a Z[i]— lattice (compare with Subsection |2.4| above). 
A Z[z]— lattice in C"^"*"^ is the free Z[i]— module generated by a basis of C™'^^ 
as complex vector space. For instance, a Z[i]— lattice is the set C 
i£^m+i q£ points in C'"^^ whose coordinates are Gauss integers. 
In these pages we consider the following two examples of lattices. First 
of all, the Z[i]— lattice 'H(^)(Z[z]) of all sequences of homogeneous polynomi- 
als F := (/i, . . . , /„) G 'H(rf) whose coefficients are Gauss integers (namely, 
fi £ Z[i][XQ, . . . , Xn]). For every complex point /? G C, we also consider the 
Z[i]— lattice 'H((i)(Z[i])(p) given by the following identity : 

W(rf)(ZH)(p) ■.= {pF : FGW(,)(Z[i])}. 

Definition 30 A non-zero point pF := (pzQ, . . . , pzjsf) G 7i(^^-^{Z[i]){p) is 
said to be C— visible from the origin if and only if the following holds : 

9cdz[j]{F) := gcdz[i\{zo, . . . , zn} G Z[z]*, 

where gcdx\i]{zo, • • • , ^^a^} is the greatest common divisor of zq, . . . ,zn in Z[i]. 

Observe that for p G M \ {0}, the class of C-visible points of 'H(rf)(Z[i])(p) is 
not equal to the class of visible points of Z^^"^^(p) (in the sense of Subsection 
mi above). 

Let P('H(^)) be the complex projective space of dimension defined above 
and let P('H(^)(Q[i])) be the A— dimensional projective space defined by the 
'^(d)(Q[i])- Let TTc : 'H(d) \{0} — > P('H(rf)) be the canonical projection. Let 
us also denote by ttc the canonical projection 

vrc : W(,)(Q[i]) \ {0} P(W(rf)(Q[i])). 

In order to introduce the notion of unitarily invariant height of a projective 



point in P('H(^)(Q[i])) (see also Subsection 2A for a comparison) we resume 



here in a very concise form the language and notation used for absolute 



values over number fields. For an introduction refer to e.g. |17, Chapter 1], 



whereas a more complete exposition of the theory of absolute values can be 



found in Artin's Algebraic Numbers and Algebraic Functions or |18]. 
Let \ ■ \u ■ Q.[i] — * ^+ be an absolute value defined on the number field 
By Q[i\u we denote the completion of Q[i] with respect to this absolute 
value I • 1,^ and by we denote the algebraic closure of Q[i];/. Finally, we 
denote by Uu the degree of over the completion of Q with respect to 
the absolute value | • |i/ : Q — > M+. 

Let MQ[j] be a proper set of absolute values of Q[i] in the sense of p7| ]. 
We assume that MQ)[i] has been chosen such that it satisfies Weil's product 
formula with multiplicities n^, i.e. for all z G Q[i] \ {0} the following holds 

n i^i^^'=i (22) 
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Let II • II : C"^"*"^ — > the canonical Hermitian norm in C™^^. Let 
S C M(Q[j] be the class of all sub-indices v £ Mq[{^ such that the absolute 
value \ ■ \i, ■ — > is a non-archimedean absolute value. 

Definition 31 For every point -kc{F) := (zq : . . . : zn) G P(W(d)(Q[i])), 
where zq, . . . , z^ £ we define unitarily invariant height of ttq^F) as : 



HaMF)) := \\F\\a I n max{| 



\ 1/2 

Zo\u,...,\zn\uV'' 



Since Weil's product formula (Equation ( ^2|) above) this quantity is well- 
defined and independent of the affine point F = (zq, . . . , zn) G 'H((i)(Q[i]) \ 
{0} we have chosen. 

Remark 32 We may also define the Northcott- Schmidt height of a projec- 
tive point 7rc(-F) E P(7Y(rf)(Q[i])) as in Subsection 2^ above. Namely, we 
may introduce 



H{MF)) ■■= \\F\\ (llmax{\zo\u,. 



\ 1/2 



Let us observe that this notion of height is essentially equivalent to Weil's 

absolute height of the projective point Trc{z) as used in p^ , JJ^ or and 
the references therein. 

From the Definition of A^,^) one easily concludes : 

^^Mpl < Ha{MF)) < ^(vrc(F)). (23) 

However unitarily invariant height is better suited for our purposes that the 
original of Northcott and Schmidt. 

Unitarily invariant height and C— visible points are closely related. 

Lemma 33 With the same notations and assumptions as above, the follow- 
ing properties hold : 

i) Let pF G 'H(^)(Z[f])(/3) be a point in that 7j[i] — lattice. Then, pF is a 
C— visible point in 7i(^^-^{'Ij[i]){p) if and only if 

\\F\\a = mm{||G||A : pG G n^a){m){p) \ {0},MpF) = MpG)}. 

a) Let pF € 'H(rf)(Z[i])(p) be a C— visible point. Then, 

\\F\\a = HaMpF)). 
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iiij Let V C f{7i^(i)) O' subset and let V := tt^^{V) U{0} be the complex 
cone over V. Let B/\(0,H) be the closed ball in C"^~^^ of radius H 
centered at the origin. Then, the following holds : 

• For every 7rc(-F) G n Pm(QH) there is a C-visible point G £ 
ynW(rf)(Z[i]) such thatiTc{G) = 7rc(F) and \\G\\a = Ha{itc{F))- 

• The number of C— visible points of that lie in V Ci 
Bj\iO,H) is four times the number of projective points in V Ci 

of unitarily invariant height at most H. 

In fact, given projective point itc{F) G P('H((f)(Q[i])) and given a C~visible 
point G G such that vrc(G) = ttc{F), the four C-visible points in 

W(rf)(Z[i]) associated to ttc{F) are {G, —G,iG, —iG} C 7i(^j_-^{Z[i]). 
Given a projective point irciF) G P('H(rf)(Q[i]), we may represent this point 
by means of some C— visible point G G 7i(^d-j{7,[i]) such that 7rc(G) = 7rc(-F). 
Thus, we may define the bit length of the projective point itc{F) as the 
bit length of this representation G. Actually, for points in W(rf)(Z[i]), the 
number of tape cells required to represent G in a Turing machine essen- 
tially agrees with the logarithm of its norm, i.e. log ||G||. As in Subsection 



2.4 above, we may define the bit length of ttc{F) as the logarithm of its 



Northcott-Schmidt height. Namely 

blMF)) := log ||G|| = logF(7rc(F)). 

As observed in Remark ^ above, Northcott-Schmidt height and unitarily 
invariant height are essentially equivalent. In Fact, the following holds for 
every 7rc(F) G P(?^(d)(Q[f]) : 

blMF)) - D^a) log%) < log HaMF)) < bl{F), 
where D(^) := max{(ii, . . . , Thus we introduce the following notion : 

Definition 34 For every ttc{F) G ¥{7i^^^{Q[i])) we define its unitarily in- 
variant bit length as the logarithm of its unitarily invariant height. Namely, 

blAiMF)) :=logFA(7rc(F)). 

5 The Non— linear Case 

5.1 Davenport's discrepancy in the Non— linear Case 



In this Subsection we apply Davenport's strategy of Subsection above to 
the Non-linear case. Given two positive real numbers -fT, e G M, let us denote 
by Ra{£,H) the compact subset defined in Equation |2l| above. Namely, 

RA{e,H) := (7r^\^{e))[J{0}) f]BAiO,H). 
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Proposition 35 With the same notations and assumptions as in Section 
above, let e G M &e two positive real numbers. Assume H > 1. Let 
(d) := (di, . . . , dn) be a list of degrees and := max{di, . . . , d^}. 
Let N{£, H) be the number of points of the Z[i] — lattice 'H((^)(Z[i]) that belong 
to Ri^{£,H). Namely, 

N{e,H) :=tt(iiA(e,^)nW(^)(Z[i])) . 

Then, the following inequality holds : 



\N{e, H) - Vol{R^{e, H))\ < -f^ — —H 



det(A(,))2 



where is the constant defined in Identity |^ of the Introduction and 

Tn := max{8,4(Z)(,) + i)}4[^^+3(n+2)2(iv+i)] ^ 



Proof- First of all, observe that R^{e,H) is the set of all systems of ho- 
mogeneous polynomials F := (/i, . . . , /„) G '^{d) that satisfies the following 
properties : 



\\F\\a<H, 

3C G C^+i, such that /i(C) 



fniO = 0, 



From Proposition ^ above we may rewrite these three properties in the 
following terms. A polynomial F G R{£, H) if and only if the following 
properties hold : 

• \\F\\l<H^ 

• 3C G such that /i(C) = • • • = /„(C) = 0, 

• 3o" G C/(n + 1), such that the following holds : 

- a(C) = eo = (1,0,...,0), 

— 3N G 7W„(C) such that det(A^) = and the following inequality 
holds 

DaiF){eo)Diagid;^^^) - \\F\\aN 

^ M2, ^ Z2 and, accordingly, let us 



From now on, let us identify 
identify 
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We also denote by U{n + 1) the corresponding real algebraic variety in 
Observe that U{n + 1) C M^" is given by a finite number of quadratic 
equations (i.e. homogeneous equations of degree 2) that may be rewritten 
as era* = Idn+i ■ 

Observe that for every a G U{n + 1), the inverse of a is the transpose 
conjugate of a (i.e. := ct*). In particular, given a := (uij) G C/(n + 1) 
and given F := (/i, . . . , /„) G "H^^^), the coefficients of cr(/j) are polynomials 
in the coefficients of fi and the entries Uij of a of degree di + 1. This is 
going to be used in the sequel. 

For every m, 1 < m < 2N + 2, let T„j be the class of all subsets / C 
{1, . . . , 2N + 2} of cardinality m. For / G Im, let J C / be a subset such 
that tt(/ \ J) = 1. Let a := (oj : j £ J) e R"*~^ be an affine real point. 
Finally, let us consider the real algebraic variety : 

Vie, H, a, I, J) C ^2N+2+2in+l)+4{n+l)^+4^ 

given by the following properties : 

A point (F, C, a, N, tiMM, U) G X M2(n+1) ^ ^2{n+lf ^ ^2(n+l)2 ^ ^4 

belongs to V{£, H,a, I, J) if and only if the following properties hold : 
r) \\F\\l-H = tl 

ii) Fj = aj, Vj G J, where F := (Fi, . . . , F2N+2) G R'^^+^ . 

Hi) F{Q = (0, . . . , 0), where F is decomposed in the list of real and imag- 
inary parts of the coefficients of the polynomials fi G . , 1 < i < n 
such that F := (/i, ...,/„), C := (Ci, ••• , C2(n+i)) G M2("+i) is iden- 
tified with the corresponding complex point C, G C"''"-'^ and F{Q is 
understood as the real and imaginary parts of the complex numbers 
F{0 G C". 

iv) a{Q = (1,0, ... ,0) G R^^""^-*^). Taking real and imaginary parts again 
these are 2(n + 1) quadratic equations. 

v) ua* = Idn+i, taking real and imaginary parts these are also 2(n + 1) 
quadratic equations. 

vi) \\F\W = f| and t^ — t^ = 0. In other words, we rewrite as quadratic 
polynomial equations the fact ts is the norm of F with respect to 
< ■,■ >A, i-e. ts = \\F\\a. 



vii) and, finally 



Da{F){eo)Diag{di - tsN - t^£^ = tl (24) 

F 



2 
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Observe that Equation ( p4D is a polynomial equation of degree at most 
max{4, 2{D^d) + 1)}, where D^a) '■= max{di, . . . , 

Now, we are in conditions to apply the upper bounds for the sum of the Betti 
numbers of a real algebraic set, simultaneous and independently obtained 
by J. Milnor (cf. R. Thorn (cf. ||), O. Oleinik (cf. |2|) and A. 

Petrovsky (cf. |^^). These upper bounds imply the following inequality : 

/?o(F(e,i?,a,/, J)) < (max{8,4(L>(d) + 1)}) ' ' ' ' , 

where /?o(^(e, H, a, /, J)) is the number of connected components of the real 
algebraic set V{e, H,a, I, J). 

Let n : M2iV+2+2(n+i)+4(n+i)2+4 — , ^2N+2 ^^le projection given by : 

n(F, C, a, N, h,t2,h, U) := F e R^''+^. 
Next, let : M^^^^ — > be the projection given by 

$/(F) := {Fi : i £ I) £ R'". 
Let rj{a) C M"^ the real line given by the conditions : 

rj{a) := {{xi : i G /) G : xj = aj,yj £ J}. 
Then, the following holds : 

^i{U{V{e, H, a, I, J))) = $/(i?A(e, H)) f| rj{a). 
As and IT are continuous mappings, we conclude 
f5o{<^i(RA{e,H)f]rj{a))) < (max{8, 4p(,) + 
With the same notations as in Subsection p.3| we conclude 

h(RA{e, H),m) < (max{8, 4p(^) + i)})27V+6(n+i)^+6 _ ^25) 
On the other hand, for every / G X^, the following holds 

where (0, H) is the closed ball in R™ of radius H, centered at the origin 
with the norm induced by || • ||a on C r2^+2. Hence, we conclude from 
the Definition of A(^) that the following holds : 

Vol{<^i(RAie,H)))<H-' ^™ 



det(A(d))2 
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With the notations of Subsection 2.3 the volume estimates V {R^{£^H),'m) 
satisfy the fohowing estimate : 

_ _ TTm(2N+2\jy' 



The following inequality is a consequence of Theorem 12 



2N+1 

\N{e,H)-Vol(RA{e,H)\ < ^ h(RA{e, H)f^+^'^V(RA{£, H), 



As H > 1, from Inequalities (|25D and (|26D we conclude : 

|iV(..i^)-F„,(fi,(..i,)|<-j^^/f-«, 

where and ©(2^+^) are the constants stated above. ■ 

The following Proposition follows by the same argument. 

Proposition 36 With the same notations and assumptions as in Proposi- 
tion 3t above, let H he a positive real number, H >1, and let N{1,H) 
be the number of systems of homogeneous polynomials in that be- 
long to the ball of radius H centered at the origin in 'Hi^d) with respect to the 
unitarily invariant norm \\ ■ Namely, 

Nil, H) := (Sa(0, H) fl (Z[i])) . 

Then, the following holds : 

_ «(2Af+2) 
|A.(l,/f)-y„,(B.(0,if))|<;jjjj^if-«. 

In particular, we conclude the following inequalities 

J^2N+2_tt2N+2 6(^^+^) u2jV+l ^ tt\ 



5.2 On C— visible ill— conditioned systems of polynomial equa- 
tions 



In this Subsection we use the notion of C— visible points (Subsection 4.4 
above) to count points in P('H(^)(Q[z])). 
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Proposition 37 Let H,e G M. be two positive real numbers, H > 1. Let 
77(8,11) be the number of points inP('H((^)(Q[i])) of unitarily invariant height 
at most H which belong to Then, the following inequality holds : 

_ Vol{R^{e,H)) ^ ^ H 

' ^ ' ^ 4C(2iV + 2) ' - ^ ' ^ 4 ' 

where _ 

^ 4det(A(,))2C(2iV + l)+ 2 

and Tn is the constant defined in Proposition above. 

In particular, the following inequality is a consequence of Corollary W^: 



■^(''^)-4C(2iV + 2)^ +L(H,iV)i^ +j, 
where is i/ie constant introduced in Corollary and 

4det(A(,))2C(2iV + l)+ 2 ' 



Proof- Assume that e > is fixed throughout this proof. We also identify 
Tii^d) = C^+^ and W(d)(Z[i]) ^ Z[i]^+^. Let F(e) C C^+i be the compact 
subset given by the following identity : 

V[e) := i?A(0, 1) fl (7rc^(S(e)) U {0}) = R^{e, 1). 

For every positive real number t £ M-|_, the following identity holds : 

tV{e) := {tF:FG V{e)} ^ R{e, t). 

In particular, for every positive real number t € M+ the following identity 
also holds : 

t'^^+^Vol{V{e)) = Vol(R{e, t)). (27) 



As in Subsection iA, for every positive real number p G M+, let 



be the Z[i]— lattice in C'^''"^ given by the following identity : 

Z[i]^+\p) := {pX :X eZ[i]^+^}. 
Let us define the finite set Ap C B/\{0, 1) given by the following equality 

Ap:=V{e)f]Z[i]'^+\p)\{0}. 

Let g{p) be the number of points in Ap. Observe that the following equality 
holds : 



g{p) := UAp) = tt (RA{e,p-')r\m''^') " 1 = N{e,p-') - 1. (28) 
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For every A G let us define the set Q{X,p) in the fohowing terms : 

n{X,p) := {XpXeVie) : gcd^[^{X) G 

where gcdx[i\ {X) is the greatest common divisor of the coordinates of X and 
Z[i]* := {1, —i} (cf. Subsection iA above for more details). 
For every m G N, m > 1, let us define : 

n{m,p):= U ^(^'/o)- 

|A|=v^ 

The following holds : 

Fact .— The following is a decomposition of Ap as a disjoint union of subsets : 

Ap=\J n{m,p). (29) 

m>l 

As n,{m, p) C Ap, this Fact is an immediate consequence of the following 
two claims : 

i) Claim 1.- Given m, m' G N, m > 1, m' > 1, if i}{m, p) n n{m' , p) ^ 0, 
then m = m' , 

ii) Claim 2.- Ap C U^>^ n{m, p). 

In order to prove Claim 1, let m,m' G N be such that m > 1, m' > 1 and 
Q{m, p) n Q{m' , p) 7^ 0. Assume 1" G Q{m,p) f] Q{m',p). Then, there are 
A, A' G such that the following holds : 

Y := XpXo = X'pX'^, 

where |A| = y/rri, \X'\ = \fm' and 



From Subsection [4.4| , these last properties mean that Xq and are C— visible 
from the origin. In particular, 

II-'^oIIa = -f^A(-'^o), II-'^oIIa = H/\{Xq). 

Now, observe that vrc(l^) = ■K£,{Xo) = 7rc(Xo)- Thus, from Weil's product 
formula we conclude the following identity : 

II^oIIa = Ha{Xo) = Ha{Y) = Ha{X'o) = WX'oWa. 

Finally, the following holds : 

||y||A := |A|p||Xo||a = |A'|p||X^||a ^v^=\x\ = \X'\ = V^, 
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and Claim 1 follows. As for Claim 2, let pX G be a point in the lattice 
'L[iY^~^^{p). With the same notations as in Subsection 4^, let A G the 
Gauss integer given by the following identity : 

A := gcdz[i]{X). 

Then, there is Xq £ such that c/cd^^f,] (Xq) G Z\i]* and X = XXq. 

In particular, 

pX = XpXo e n{x,p) c ni\x\,p). 

Taking m := |Ap > 1, Claim 2 follows. 

For every integer number m G N, m > 1, let f{mp) be the number of points 
in 0,{m,p). Observe that this function is well-defined for fixed p. Then, 
Equality ( p9| ) above, yields the following equality : 

oo 

9{p) ■■= f{mp). 

Applying Mobius inversion formula (see Theorem 270 in |14], for instance) 
it follows that 



m=l 



where p is Mobius function. From Theorem 287 of |ll4|| (cf. Equation (|ll| 
above) we have 



2N+2 fi N Vol{V{e)) -sr^ / 2N+2 ^ ^ ^ ^ JjW)\ 

^ \ ' m=l ^ ' 



Hence, we conclude 



00 



^ ' 171=1 

(30) 

As in the Proof of Proposition 18, let Np be the set of integer numbers given 
by the following identity 

Np := {m G N : I < m, mp < 1}. 

Let us denote by Sp the sum of the terms of the left hand series of Equation 
(30) whose indices are in Np. Namely, 

^p--= E ^^\impr'-'9{mp)-Vol{V{e))\ 
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As in the Proof of Proposition ^ above, assume m E Np. Then, replacing 
{mp)~^ by we conclude from Identities ( p7|) and ( p8|) above the following 
identity : 

\{mpf''+^g{mp) - Vol{V{e))\ = -^^|iV(e,F) - {I + Vol{R^{e,H))\ 

(31) 



From Proposition 35 above, we have : 



1 _ _ r/v6(2^+2) 

:\N{e,H)-Vol{R^{e,H)\ < ^ 



Replacing back H by {mp)~^, Equation (31) becomes the following inequal- 
ity : 



limpf^'+'gimp) - Vol{V{e))\ < (mp)^^^^-^ + {mpf''+\ (32) 
We conclude : 



meNp 

Then, 



^ E I det(A(,))^ + ^("^^^ ) • ^^'^ 



''-C(2iV + l)det(A(,))2^J- ^ 



and we conclude 



^^-C(2iV + l)det(A(,))^+^ • (''^ 
As in the Proof of Proposition let be the class given by 

N'^ := {m G N : 1 < m, mp > 1}. 

Observe that if m G N^, then g{mp) = 0. 

Hence, let S'p be the sum of the terms in the left hand series of Equation 
whose indices are in N^. Namely, 

meN' 



Then, we have 

^ E ^T2\^^l(yie))\ < pVol{V{e)) (35) 
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Now, we combine Equations ( p^ ) and ( pq ) with Equation (|3^) to conclude 



Finally, replacing p ^ by -fT in this Equation and using the identity described 
in Equation ^ we conclude : 

1 Vol(RA{e,H)) T^6(2A^+2) ^^^^ yo/(fiA(g, ^ )) 

'^^ ^ C(2iV + 2) ' - C(2A^ + l)det(A(rf))2 ^ ^ H 

Finally, from Lemma ^ above we have ^^-^^ — - = AA(e, H) and Proposition 
37 is achieved. ■ 

The same arguments yield the following Proposition 

Proposition 38 Let H,e £ M be two positive real numbers, H > 1. Let 
A/'(l,i?) be the number of points in P('H(^)(Q[i])) of unitarily invariant 
height at most H . Then, the following inequality holds : 

^) - 4C(2iV + 2) ' - + T' 

where 

e(27V+2) 



4det(A(d))2C(2Af + l) 2det(A(rf))' 
In particular, the following inequality also holds : 

K2N+2 



4det(A(rf))C(2iV + 2) 



5.3 Proof of Theorem g 

We are now in conditions to prove Theorem Q as stated at the Introduction. 
Proof of Theorem 0.- First of all, observe that the probability that a random 
choice of a system of polynomial equations F S P(7^(rf)Q[i]) of unitarily 
invariant height at most H satisfies iinorm{F) ^ 1/e is given by the following 
quotient : 

Jf{e,H) 

Let A be the quantity : 



A :-- 



4det(A(rf))C(2iV + 2) 
From Propositions pq and conclude : 
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H 



7^(1, H) > AH 



2N+2 



L{l,N)H 



2N+1 



+ 



H 



where := n^{n+l)N[N —1)'D(^^-^ is the constant introduced in Theorem 

^ of the Introduction. Then, we conclude 



< m)r + 



H-g{N, H)' 



where 



and 



A 



T{e,N,H) := £^^[{d)]g{N, H) + 
Simphfying these expressions we conclude : 

q{2N+2) 

G{N,H) = C{2N + 2) 



L{e,N) 



+ 



1 



+ 



det(A(rf)) 



det(A(rf))C(27V+l)i^2iV+2 i^27V+2^27V 

On the other hand, the following identity holds : 



+ 2 



A 



C(2iV + 2) 



rA,6(2iv+2) 



Hence, we conclude 

J^{£,N,H) =C{2N + 2 



det(A(rf))C(2iV + l)i^27V+2 



+ 2s^^[{d)] 



det(A(rf)C(2iV + l)^2JV+2 



+ 



+C(2iV + 2 

As det(A(^)) < 1 and H > 1, we may conclude : 



where 



Mil,H) 
g{N): 



C{2N + 2)-'^H -g{Ny 

q{2N+2) I 



+ 



det(A(rf))C(2iV + l)K2N+2 K2N+2 



+ 2, 



and 



>■ det(A(,)C(2iV+l)K2;v+2 i^27V+2 
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6 Applications : On the average precision required 
to write down an Approximate Zero 



In this Section we discuss some applications of the estimates of Theorem ^ 
to compute the average precision required to write down approximate zeros 
for the affine Newton operator. We start by recahing some elementary facts 
about these ideas. 



6.1 Some well— known facts about AfRne Approximate Zero 
Theory. 

Let (d) := {di, . . . , dn) be a list of degrees, and let P(^) be the set of all 
sequences F := (/i,...,/„) of multivariate polynomials that satisfy the 
following properties for every i, 1 < i < n : 

i) /i gC[Xi,...,X„] and 

ii) deg(/i) < di, 

We may easily identify V^^^-^ and Ti-(d) by the obvious homogenization opera- 
tor. For every F £ 'P{d)j we also denote by F its homogenization F £ 'H[d)- 
In this Subsection we study the affine counterpart of multivariate Newton's 
operator. Let F := (/i, . . . , /„) S Vi^d) be a sequence of multivariate polyno- 
mials. Let V{F) C C" be the set of common zeros into the complex affine 
space C*, namely 

y(F):={CGC'^ : /i(C) = • • • = /„(C) = 0}. 

Let C, € V{F) be a smooth zero of system F, namely the jacobian matrix of 
-F at C is non-singular, i.e. DF{C,) G GL{n,C), where 

\ J / l<t,j<n 

Let Np be the multivariate Newton operator defined by F. Namely, 

^1 

NFiXi,...,Xn):= \ l \ -DF{Xi,...,Xn)-' 



h{Xi,... 


) Xn 


fn{Xl, . . . 


, Xn 


p, n 





, Xn 

Next, we recall Smale's Approximate Zero Theory (cf. |3^, |2^], |2^], p9| ], 
[30 1, 1 32 1, 1 31] and the compiled version in [^) as foundation of Numerical 
Analysis. 

With the same notations and assumptions as above, an approximate zero 
of the system F with associated zero is a point z E C" such that the 
sequence of iterates of the Newton operator is well-defined and converges 
quadratically to C- Namely, an approximate zero of the system F with 
associated zero C £ ^(-^) is an affine point z £ C" such that the following 
properties hold : 
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• The jacobian matrix of F at z is a regular matrix. Namely, DF{z) G 
GL(n,C). 

• The following sequence is well-defined : 

zi := Nf{z) G C", and Zk := Npizk-i) for k>2. 

• For every A; G N, k > 1, the following inequality holds : 

\\zk-C\\<^\\z-C\\ 

With the same notations as above, we define the quantity 7 in the following 
terms : 



l{F, C) := sup 

fc>2 



1 

fc-1 



where the norm on the right hand side of this identity is the norm of the 
multilinear operator 

This quantity yields a locally sufficient condition for having an approximate 
zero. 



Theorem 39 (7— Theorem, [ p3|| ) With the same notations as above, Let 
F G 'P(rf) be a regular sequence of polynomials that defines a smooth, zero- 
dimensional affine algebraic variety V{F), and let C G C" a smooth zero of 
system F. Let z G C" be an affine point satisfying the following inequality : 

||C-^ll7(/,C)<^^. 

Then, z is an approximate zero of the system F with associated zero C,. 

For every affine point C := (Cii • • • > Cn) S C^, let us denote by C £ IPn(C) the 
corresponding projective point given by the following identity : 

C:= (l:Ci:---:Cn)GlPn(C). 

Combining Lemma 7 and Theorem 2 of Chapter 14 in for every Q G C", 
such that -F(C) = and DF[C,) G GL{n,C), the following inequality holds : 
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6.2 On the number of approximate zeros of given precision 

Let us assume now the computational hypothesis (cf. also Q) : Assume 
your are given a system of polynomial equations F S V(^ij_) ^^'^ you want 
to compute an approximate zero z £ C" of system F with associated zero 
E C". From a computational point of view, you look for approximate zeros 
z that could be written by a computer (or equivalently that can be written 
as a list of symbols in a finite alphabet). In standard numerical analysis 
procedures, the computational hypothesis becomes now the assumption that 
z G Qii]'^. 

Given a point z G Q[i]", we define the precision of z as the bit length of its 
denominator. Namely, let z € Q[i]"' be a point whose coordinates are Gauss 
rationals. Then, there are g G Z \ {0} and (zi, . . . , z„) G such that 



z :- 



£i 



and q is of minimal absolute value satisfying this property. The precision 
Pr{z) of z G Qfi]" is defined to be Pr{z) := max{0, log2 g}. 
The following statement gives upper and lower bounds for the number of 
approximate zeros of given precision that satisfy the 7— Theorem above. In 
other words, we show upper and lower estimates for the behaviour of the 
function : 

NmiC) ■•= ^{z G Q[i]" : \\z - Cll < l^^y^' A Pr{z) = log m}. 
Observe that Nm{C) equals the cardinal of the set : 

7(F,C) ' ' 



Proposition 40 With the same notations as above, the following properties 
hold : 



i) For every integer number H such that 



H < 



the following holds 



,3-^7. 



H 



m=l 
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ii) For every m G N such that 



Hi := — <m<H2:-- 



3-V7J - - 3-^/7' 
the following holds 

< N^iC) < 1. 

In particular, 

m=Hi ' 

Hi) For every m E N the following inequalities hold : 

f m(3-V7) ^\ ^" _ f m(3-V7) ^/2^\ ^" 



Proof- Item ui) follows from an elementary argument as those of [^] or 
the more classical Blichfeldt-Minkowski estimates. Similar estimates could 
also be obtained by means of Davenport's Theorem (Theorem [l^ above). 
We left the proof of item iii) for the reader. As for items i) and ii) we make 
use of the following estimate, which is a translation of the "Gap Principle" 



of |21] to our context (cf. also p5|]). 

Given z,z' G Q[i]"' two different approximate zeros of system F with asso- 
ciate zero C that satisfy the 7— Theorem, and such that mz, mz' G the 
following chain of inequalities hold 

1 ,,mz mz',, ,, ^„ „ , 3 — \/7 

- < < \\Z - C + \\Z - C < -7^- 

m m m ^[r^Q) 

Hence, we conclude m > -^"-^ and the following holds for every m such 
that 1 < m < : 

NmiO < 1- (37) 

Item i) follows from the following Remark. Assume that H satisfies the 
following inequality 

Let z, z' G Q[i]" be two points such that there are d, d' £ N \ {0} satisfying 
the following properties : 

• The coordinates of the points dz, d' z' £ Z[i]" are Gauss integers. 
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• z^z' are approximate zeros of system F with associate zero C that 
satisfy the 7— Theorem above. Namely, assume the following holds 

• d,d' < H. 

Then, the following chain of inequalities hold : 

2,<||._.'||<||._cll + l|.'-CII<^. 

In particular, we would conclude dd! > ■ In other words, '\i d < H, 

3— V 7 

then d! > H and item i) follows. Item ii) follows from the previous Remark 
and Inequality ^ above. ■ 



Item in) of Proposition 40 immediately yields the following 

Corollary 41 With the same notations and assumptions as above, for every 
p G N such that the following inequality holds : 



P > log7(F,C) +log 



+ 1, 



there are approximate zeros z £ Q[i]" of system F of precision p with asso- 
ciated zero (. 



The following statement follows from the previous Corollary, Identity 36 and 
our estimates on the probability distribution of finorm- 

Corollary 42 With the same notations as above, there is a universal con- 
stant C3 > ( C3 < 20 j such that the following holds : 

Let (d) := {di, . . . , dn) be a list of degrees and let h,w £ M. be two positive real 
numbers. Assume that h is big enough to satisfy the following inequality : 

h > C3A^2(logAf + logD(d)) +logu;. 

Then, for every system of polynomial equations F S P('P((^)(Q[z])) of bit 
length at most h, the required precision Pr{F) to write down an approximate 
zero of F satisfies 

Pr{F) < 0{n log2 + log^ w), (38) 
with probability at least 

2 

1 . 

w 
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